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ABSTRACT 


Empirical  comparison  of  the  two  segment  ma;:imum 
likelihood  f-k  spectra  with  the  fast  frequency  domain 
beam  f-k  spectra  (FKCOMB,  FKPLOT)  shows  that  the  latter 
is  more  suitable  for  the  separation  of  multiple  signals 
if  the  amplitudes  of  the  signals  differ  considerably. 
This  advantage  of  the  beaming  process  is  attributed 
primarily  to  the  stripping  procedure,  which  has  not 
been  developed  for  the  maximum  likelihood  spectra. 

The  maximum  likelihood  f-k  spectra,  on  the  other  hand, 
are  less  sensitive  to  the  array  response  and 
easier  to  interpret,  especially  if  the  array  used 
cortains  only  a few  elements,  since  the  sidelobes  of 
the  array  response  are  more  confusing  on  the  FKCOMB 
output  of  this  case.  Both  processes  require  a good 
signal  to  noise  ratio  (-2)  for  successful  application. 

Composited  recordings  of  long  period  Rayleigh  waves 
recorded  at  LASA  were  utilized  for  the  comparison  using 
various  subarray  configurations  and  signal  and  noi.  e 
levels . 
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INTRODUCTION 


The  purpose  of  this  report  is  to  compare  the 
performance  of  two  ways  r T estimating  the  frequency 
wavenumber  structure.  The  first  method  is  a fast 
frequency  domain  beamforming  program  which  was  devel- 
oped at  the  Seismic  Array  Analysis  Center  by  E.  Smart. 
Of  the  various  programmed  versions  of  this  method,  we 
used  FKPLOT , which  forms  the  frequency  wavenumber 
spectra  at  a selected  frequency.  A later  version 
(FKCOMB)  searches  for  the  frequency  with  the  highest 
signal- to-noise  ratio  and  automatically  determines  the 
azimuths  and  velocities  of  the  prominent  peaks  in  the 
f-k  plane.  Detailed  descriptions  of  these  programs 
and  the  method  can  be  found  in  SAAC  Report  No.  9; 
and  in  papers  by  Smart  and  Flinn  (1971)  and  Mack  and 
Smart  (1971).  An  important  feature  of  this  method  is 
a procedure  called  "stripping"  which  subtracts  the 
most  prominent  spectral  peaks  and  the  associated  side 
lobes  from  the  f-k  spectra  to  reveal  any  hidden  weak 
signals.  Blandford,  Cohen,  and  Woods  (1973)  have 
shown  that  stripping  is  the  first  iteration  of  an 
approximation  to  a two-signal  maximum  likelihood  pro- 
cessor. 

The  second  method  used  is  the  maximum  likelihood 
f-k  analysis  described  by  Capon  (1969),  Lintz  (1968) 
and  McCowan  and  Lintz  (1968).  The  output  obtained  by 
this  method  is  the  power  passed  by  a set  of  maximum 
likelihood  filters  which  were  designed  to  pass  a plane 
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wave  corresponding  to  a particular  point  in  the  f-k 
plane  and  reject  the  rest  of  the  waves  present.  The 
computation  of  the  maximum  likelihood  spectra  uses  an 
estimate  of  the  power  spectral  matrix  of  all  the  sen- 
sors. In  the  method  originally  proposed  by  Capon  et 
al.  (1967),  this  power  spectral  matrix  was  estimated 
by  the  Fourier  transformation  of  many  time  segments  of 
data,  forming  all  the  auto-  and  cross-products  for 
each  segment,  and  finally  averaging  over  many  time 
segments.  This  procedure  gives  the  necessary  statis- 
tical stability  to  the  estimates.  The  procedure 
described  can  be  used  for  long  noise  samples,  but 
difficulties  arise  when  dealing  with  earthquake  signals. 
The  duration  of  these  signals  is  relatively  short  and 
the  stability  of  the  spectral  estimates  has  to  be 
sacrificed.  For  the  computations  given  in  this  report 
only  two  segments  of  128  seconds  duration  were  used; 
we  call  the  output  the  "two  segment  maximum  likelihood 
f-k  spectra"  (TSML  f-k  spectrum) . 

In  comparing  the  two  methods  we  examined  the  fol- 
lowing features: 

1.  Ability  of  the  two  processes  to  separate  two 
events  close  to  azimuth  of  arrival  and  nearly  equal 

in  amplitude,  with  high  S/N  ratio. 

2.  Ability  to  separate  two  events  close  in 
azimuth  and  unequal  amplitude,  with  high  S/N  ratio. 

3.  Sensitivity  of  the  two  processes  to  noise. 

4.  Effect  of  reducing  the  array  size  and  the 
number  of  sensors. 
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Theoretical  analyses  of  the  two  methods  have  been 
made  for  idealized  cases  which  show  that  the  maximum 
likelihood  method  is  in  many  ways  superior  to  the  fast 
beamed  f-k  spectra.  The  resolution  should  be  consid- 
erably better,  the  TSML  spectra  in  the  noiseless  case 
should  be  independent  of  the  array  response,  and  two 
uncorrelated  signals  at  different  azimuth  should  be 
completely  separable  (Lintz  and  Woods,  1972).  It  is 
not  clear,  however,  how  these  processors  would  perform 
in  the  presence  of  noise,  and  what  effect  accidental 
correlations  of  the  two  signals  within  the  short  time 
windows  would  have.  In  the  following  we  will  define 
these  events  to  be  "separated"  if  two  peaks  appear  in 
the  TSML  output  which  correspond  to  azimuths  and 
directions  of  the  known  signals  and  if  these  peaks  are 
at  least  2-3  dB  above  the  general  background.  The 
same  definition  is  used  for  the  FKPLOT  (original  and 
stripped)  outputs.  An  additional  criterion  requiring 
that  the  peaks  in  the  f-k  plane  should  correspond  to 
realistic  Rayleigh  wave  velocities  will  also  be  used 
occasionally.  The  velocity  3.6  km/sec  corresponding 
to  average  phase  velocity  of  Rayleigh  waves  around 
30  sec  period  is  marked  by  a circle  on  all  f-k  plots. 
Infinite  velocity  is  marked  by  a cross  at  the  center 
of  the  plot. 

The  contouring  used  in  the  presentation  conforms 
to  that  of  the  FKPLOT  output,  that  is,  the  first  con- 
tour is  at  1 dB  below  the  peak,  the  second  3 dB  below 
peak,  and  the  following  ones  are  separated  by  3 dB. 


-3- 


i 


The  TSML  outputs  are  similarly  contoured,  although  the 
program  gives  the  output  at  1 dB  intervals.  The  repro- 
ductions of  the  printer  plot  shows  all  these  subinter- 
vals and  they  are  occasionally  indicated  by  dashed 
contours.  Velocities,  azimuths,  S/N  ratios  and  F- 
sfatistics  computed  by  FKPL01  are  shown  in  Table  II. 

In  testing  the  two  processes  we  did  not  compare 
their  computational  efficiency,  since  while  the  pro- 
duction-oriented FKPLOT  was  considerably  optimized 
with  respect  to  running  time,  the  research  programs 
made  available  to  us  by  Texas  Instruments  for  computing 
the  TSML  spectra  were  not.  Therefore,  it  would  have 
been  unfair  to  compare  the  running  times  of  the  two 
processes,  although  the  simplicity  of  the  FKPLOT 
beaming  piucess  makes  it  likely  that  even  with  the 
inclusion  of  stripping  it  would  be  faster  than  com- 
puting the  maximum  likelihood  spectra. 
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DATA 


Throughout  this  report,  long -period  vertical 
component  data  for  LASA  have  been  used.  Because  of 
data  quality  problems,  some  sensors  had  to  be  omitted. 
For  parallel  runs  of  the  two  processors  the  number  of 
sensors  used  was  kept  the  same  to  make  the  comparisons 
valid.  The  elements  Cl  and  C2  were  omitted  at  all 
times  because  of  gain  and  data  quality  problems,  and 
element  D2  was  also  omitted  when  Event  5 was  present. 

To  test  the  processors  six  events  were  selected.  The 
epicentral  data  fox  these  events  are  shown  in  lable  I. 
Three  of  the  events  occurred  in  the  western  Pacific 
region  with  Rayleigh  waves  arriving  at  LASA  from  the 
WNW  direction.  The  three  other  events  are  roughly 
south  of  LASA.  The  Rayleigh  wave  train  for  the  Baja 
event  (.Event  4)  is  short  and  fits  entirely  within  the 
256  second  window  used  in  the  calculations.  In  order 
to  make  the  results  mutually  comparable,  the  f-k  spec- 
tra were  computed  for  a common  frequency  fQ=. 04688  cps. 
(T=21 . 3) . In  order  to  eliminate  undesirable  leakage 
of  energy  due  to  short  time  windows,  the  windows  were 
selected  by  visually  estimating  the  dominant  frequency 
along  the  wavetrain  and  placing  the  window  at  the 
portion  of  the  wavetrain  where  the  chosen  reference 
frequency  fQ  was  dominant.  Subsequent  exact  Fourier 
analysis  showed  that  the  dominant  frequency  was  suf- 
ficiently close  to  f to  rule  out  excessive  leakage 
from  other  frequencies  which  would  result  in  false 
velocities  (Smart,  1971).  The  length  of  time  windows 
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used  in  FKPlOT  was  2 56  seconds.  The  TSML  method  uses: 
the  same  length  of  record  divided  into  two  time  seg- 
ments. figures  1 to  6 show  plots  of  all  traces  of  the 
events  used  with  the  time  windows  marked.  All  western 
Pacific  events  show  considerable  variation  in  the 
surface  wave  forms  across  the  array,  which  indicates 
multipathing.  This  is  somewhat  less  pronounced  in 
New  Britain  event  (Event  2).  The  Baia  California 
event  (Event  4)  lias  a short  surface  wavetrain.  The 
Argentine  event  (Event  5)  and  the  Guatemala  event 
(Event  6)  show  considerable  multipathing,  as  evidence., 
by  beats  on  the  wavetrain.  All  events  chosen  were  oi 
a sufficiently  high  magnitude  for  the  signal  to  be 
considerably  above  the  background  noise  level.  Static 
tical  measures  of  S/N  ratio  computed  during  the  FKPL  1 
runs  are  measures  of  the  deviation  of  the  input  from  a 
single  plane  wave  rather  than  that  of  ratios  of  signal 
level  to  background  noise. 

Figure  7 shows  the  power  spectra  of  the  events 
within  the  signal  windows  chosen,  averaged  over  all 
sensors  and  the  two  time  segments  used  for  TSML  ana- 
lysis. The  chosen  reference  frequency  fQ  is  at  or 
very  near  to  the  spectral  peak  for  all  events.  Small 
differences  should  not  affect  the  results  presented. 
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F-K  ANALYSES  OF  DATA  USING  THE  FULL  ARRAY 


General 

In  the  following,  the  analyses  of  various  signals 
and  signal  and  noise  combinations  is  given,  using  all 
the  available  vertical  elements  of  the  LASA  long  period 
array.  The  results  are  representative  of  the  capa- 
bilities of  a large  array,  which  is  not  like.y  to  be 
duplicated  in  the  near  future.  In  the  latter  part  of 
the  report  only  parts  of  LASA  will  be  used,  reducing 
the  resolution  but  at  the  same  time  giving  results 
more  indicative  of  the  performance  of  smaller  arrays 
presently  being  contemplated. 

Analyses  of  individual  signals  and  noise  samples 

Figure  8 shows  the  f-k  spectra  obtained  for 
Event  1,  the  FKPLOT  output  with  the  stripped  region  in 
the  center,  the  TSML  f-k  spectrum  on  the  right.  The 
TSML  output  shows  a peak  at  azimuth  347,  with  velocity 
s-  3.6  km/sec.  The  FKPLOT  output  shows  a secondary 
peak  at  a slightly  different  azimuth  on  the  stripped 
f-k  plot,  which  is  probably  due  to  a multipath  arrival. 
In  general  the  stripped  peaks  will  be  rejected  for 
false  events  because  they  will  have  low  values  for  the 
"F"  statistic. 

Figure  9 shows  f-k  spectra  for  Event  2.  The  TSML 
output  shows  a well-defined  peak  showing  no  more  multi- 
pathing  than  Event  1.  The  FKPLOT  output  shows  a much 
broader  peak  at  roughly  the  same  azimuth,  with  high 
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S/N  ratio,  evidently  modified  by  the  array  response. 
The  stripped  output  shows  a spurious  peak  at  a low 
S/N  ratio,  which  is  probably  meaningless. 

Event  3 is  more  spread  out  in  azimuth  than  Event 
2,  as  evidenced  by  the  TSML  output  (Figure  10).  The 
difference  is  not  evident  on  the  FKPLOT  output,  since 
this  technique  has  a lower  resolution.  The  stripped 
version  again  shows  spurious  peaks  at  low  S/N  ratio. 


Event  4 shows  a single  peak  on  th 
on  Figure  11;  the  FKPLOT  outpu.  again 
peak  on  the  original  f-k  plot  and  two 
at  moderate  S/N  ratio  on  the  stripped 
these  peaks  are  multipath  arrivals  it 
on  the  TSML  output. 


e TSML  f-k  plot 
shows  a broad 
spurious  peaks 
version.  If 
is  not  evident 


Event  5 has  broader  peaks  on  both  kinds  of  f-k 
spectra  indicating  multipathing  or  scattering  of 
surface  waves  (Figure  12).  This  is  again  quite  plau- 
sible considering  the  travel  path  of  these  waves 
across  the  Andes  Mountains  and  Central  America.  The 
stripped  version  shows  a peak  at  low  S/N  ratio, 
possibly  due  to  a multipath  arrival.  Similar  conclu- 
sions can  be  reached  for  L\cnt  6 shown  in  Figure  13. 


The  above  events  are  thought  to  be  fairly  repre- 
sentative of  real  earthquake  signals  with  respect  to 
the  azimuthal  spread  of  Rayleigh  waves.  Real  earth- 
quake signals  are  not  plane  waves  and  theoretical 
evaluation  of  the  performance  of  processors  using  plane 
wave  models  can  be  misleading.  It  is  important  to  test 
the  processors  using  real  signals. 
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Figure  14  shows  the  f-k  plots  of  a noise  sample 
used  for  superposition  on  the  various  inputs  to  be 
analyzed  to  simulate  noisy  signals  Although  it  shows 
some  directionality,  it  is  thought  to  be  fairly  repre- 
sentative of  noise  in  general.  Visually  the  noise 
seems  to  be  uncorrelated  between  sensors.  The  noise 
f-k  spectra  are  included  here  for  comparison  with  the 
f-k  analyses  of  the  noisy  signals. 


ANALYSIS  OF  COMBINED  EVENTS 

In  order  to  simulate  two  signals  which  arrive 
simultaneously  at  the  array  the  signals  of  two  events 
are  added  within  the  windows  previously  used.  By 
combining  the  events  at  various  relative  amplitude 
levels  and  adding  noise  one  can  test  the  sensitivity 
of  the  two  processes  to  different  signal  amplitudes 
and  to  varying  noise  levels.  The  signal  levels  were 
scaled  to  maximum  amplitudes  occurring  within  the 
selected  windows  at  the  center  array  element  for  each 
event.  This  rough  scaling  is  sufficient  for  our 
purposes,  although  scaling  the  spectral  amplitudes  at 
the  preselected  reference  frequency  f would  be  more 

appropriate.  Various  cases  were  tested  and  the  results 
are  given  below. 

It  should  be  noted  that  the  additional  signal  will 
cause  the  statistics  which  are  indicators  of  the  S/N 
ratio  (F  statistics  and  S/N  ratio  given  by  FKTLOT)  to 
decrease.  Therefore,  these  indicators  will  lose  their 
original  meanings.  Low:  values  of  these  parameters  do 
not  indicate  the  presence  of  high  ambient  noise.  Ii 
one  event  is  removed  by  stripping,  these  indicators 
can  increase  in  value  since  tiie  energy  seems  to  come 
from  fewer  directions.  This  phenomenon  is  discussed 
by  Blandford  ( 1972)  . 
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Case  1 


The  two  signals  are  equal  in  amplitude;  no  noise 
is  added. 

The  f-k  spectra  for  ti  e combination  of  events  1 
and  2 are  shown  in  Figure  15.  The  FKPLOT  output  shows 
that  without  stripping  there  are  only  faint  indications 
of  the  presence  of  Event  1,  but  stripping  removes 
Event  2 effectively  and  Event  1 then  becomes  visible. 
The  S/N  ratio  for  both  cases  is  fair.  The  TSML  output 
records  both  signals,  with  Event  1 about  5 dB  below 
Event  2.  Thus  both  processes  resolved  the  signals 
successfully. 

The  combination  of  Events  1 and  3 are  well  re- 
solved on  both  the  TSML  and  the  unstripped  FKPLOT 
outputs.  On  the  TSML  output  Event  is  4 dB  above 
Event  3 while  the  FKPLOT  peaks  are  on  the  same  level. 
Stripping  leaves  Event  3,  vhe  weaker  signal.  Again 
both  processes  performed  eoually  well  (Figure  16). 

Figure  17  shows  that  Events  t and  5 are  not  well 
separated  on  the  TSML  output,  the  second  e^ent  is 
indicated  by  a bulge  in  the  -2  dB  contour.  Only  this 
secondary  dashed  contour  separates  the  two  events. 

The  FKPLOT  outp  >t  shows  only  a large  maximum  including 
both  events.  After  stripping  Event  5 dominates.  The 
indicated  S/N  ratio  is  small  due  to  the  large  azimuthal 
spread.  Both  processes  performed  equally  well,  re- 
sulting in  a doubtful  separation  if  no  other  informa- 
tion is  available. 
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Events  4 and  6,  however,  are  well  separated 
(Figure  18)  by  both  processes.  Stripping  leaves 
Event  6 after  Event  4 is  removed.  The  TSML  output 
shows  two  well-defined  maxima  at  about  the  correct 
azimuths.  S/N  indicators  are  moderately  low  for  the 
original  FKPLOT  output,  somewhat  higher  on  the  strip- 
ped output,  ihe  reasons  for  this  have  been  previously 
stated  in  the  discussion  above. 

Summarizing  the  above  results  it  seems  that  both 
processes  perform  equally  well  if  the  two  events  have 
roughly  the  same  amplitude  within  the  window  and  the 
S/N  ratio  is  high. 

Case  2 

One  signal  is  scaled  to  be  three  times  the  ampli- 
tude of  the  other;  no  noise  was  added. 

3ecause  the  scaling  was  done  using  maximum  ampli- 
tudes and  not  spectra,  the  actual  amplitude  ratio  is 
somewhat  different  from  3:1. 

The  combination  of  Events  2 and  1 with  the  ampli- 
tude ratio  3:1  could  not  be  resolved  by  the  TSML 
processor  but  can  be  well  resolved  by  FKPLOT  if  strip- 
ping is  used.  The  actual  amplitude  ratio  at  the 
reference  frequency  f is  higher  than  3:1  as  shown  in 
the  previous  section  (Figure  19)  , the  discrepancy 
resulted  from  the  rough  scaling  procedure  used. 

For  the  combination  of  Events  3 and  1 (Figure  20) 
with  3:1  amplitude  ratio,  (which  is  actually  somewhat 
less  at  fQ)  both  processes  separate  the  two  events 
quite  well  if  stripping  is  used  for  FKPLOT. 
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The  superposition  of  Events  4 and  5 (Figure  21) 
shows  Event  4 with  high  S/N  on  the  unstripped  FKPLOT 
output  and  Event  5 with  moderate  S/N  on  the  stripped 
output.  TSML  fails  to  separate  the  two  events  with 
certainty,  only  a bulge  in  the  -9  dB  contour  indicating 
the  presence  of  a second  signal.  FKPLOT  in  this  case 
performed  considerably  better. 

The  combination  of  Events  4 and  6 (Figure  22) 

« 

shows  separation  of  events  for  both  processes.  The 
Event  6 peak  on  the  TSML  output  is  about  6 dB  down 
from  the  Event  4 peak,  barely  detectable  and  at  not 
exactly  the  right  azimuth  but  separated  more  from 
Event  4 than  the  original  azimuth.  After  stripping, 
FKPLOT  yields  the  Event  6 peak  with  a reasonable 
certainty  but  somewhat  farther  from  the  Event  4 azimuth 
than  the  original.  In  this  example  we  feel  that 
TSML  performed  somewhat  better,  since  the  azimuth  of 
Event  6 is  better  on  the  TSML  output  although  it  is 
not  much  above  the  background  level.  However,  the 
peak  on  the  stripped  FKPLOT  is  about  4 dB  higher  than 
the  general  background,  but  the  azimuth  is  wrong. 

Overall,  if  two  signals  are  not  on  the  same 
amplitude  level,  FKPLOT  with  stripping  may  out-perform 
the  two  segment  maximum  likelihood  method. 

Case  3 

Events  combined  with  equal  maximum  amplitudes, 
noise  added.  S/N  ratio  2. 
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For  this  case  the  noise  sample  selected  was 
added  to  the  combined  events. 

The  combination  of  Events  2 and  1 with  noise 
added  (Figure  23)  does  not  separate  well  on  the  TSML 
frequency  wavenumber  spectra  where  a peak  in  the  ENE 
direction  is  only  slightly  below  the  level  of  the  very 
weak  peak  corresponding  to  Event  1.  The  separation  is 
better  on  FKPLOT  including  the  stripped  version. 

Both  processes  fail  for  the  combination  of  Events 
3 and  1 (Figure  24) . 

Comparing  the  outputs  for  the  superposed  Events  4 
and  5 (Figure  25)  it  seems  that  both  processes  fail  to 
separate  the  two  events  with  certainty.  The  FKPLOT 
output  has  a broad  peak  containing  both  events,  after 
stripping  one  obtains  a peak  coinciding  with  the 
direction  of  Event  4,  but  within  1 dB  of  many  other 
similar  peaks.  The  TSML  output  shows  a complicated 
picture,  the  main  peak  is  in  a direction  considerably 
off  the  direction  of  Event  4,  possibly  influenced  by 
the  noise  added. 

Events  4 and  6 seem  to  separate  well  on  both 
FKPLOT  and  TSML  outputs  (Figure  26) , but  the  azimuth 
corresponding  to  Event  4 is  considerably  off  on  the 
TSML  output.  This  can  be  caused  by  a noise  peak. 

Case  4 

Events  combined  with  equal  maximum  amplitudes, 
noise  added,  S/N  ratio  1. 
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The  results  of  this  test  are  shown  in  Figures  27 
through  30.  In  general  it  can  be  said  that  both  FKPLOT 
and  the  TSML  process  break  down,  showing  many  peaks  at 
the  wrong  azimuths  and  velocities,  which  have  apparently 
no  relationship  to  the  actual  signals  present,  while 
also  revealing  occasional  signal  peaks  which  cannot  be 
distinguished  from  false  peaks.  Thus  both  f-k  processes 
require  a signal  to  noise  ratio  of  about  2 to  perform 
adequately. 


F-K  ANALYSES  WITH  THE  LASA  I'  RING  OMITTED 


General 

The  omission  of  the  F ring  reduces  the  array 
diameter  from  200  to  125  km  (Figure  72)  with  a corre- 
sponding reduction  of  resolving  power.  All  runs  were 
repeated  with  the  reduced  array,  including  the  indi- 
vidual event  and  noise  samples,  to  give  a basis  of 
reference  for  comparison  of  various  signal  and  noise 
combination  f-k  plots  using  this  array  configuration. 

Individual  signals  and  noise 

The  f-k  plots  of  individual  signals  are  shown  in 
Figures  31  through  56.  There  are  no  basic  differences 
between  the  plots  shown  and  those  computed  using  the 
full  array  except  for  the  reduction  in  the  resolution 
for  both  processes,  although  the  resolution  should 
not  reduce  as  much  in  the  case  of  TSML  according  to 
theory.  The  sidelobe  pattern  is  naturally  different 
for  this  array  configuration,  but  otherwise  the  removal 
of  F ring  did  not  make  much  difference.  Event  5 
(Figure  35)  shows  a peculiar  double  arrival  pattern, 
not  seen  before.  The  noise  sample  analysis  are 
shown  in  Figure  37. 
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ANALYSIS  OF  COMBINED  EVENTS 

Except  for  the  array  configuration,  all  cases 
are  identical  to  those  mentioned  previous-^. 

Case  1 

Events  are  equal  in  amplitude;  no  noise  added. 

As  in  the  case  of  the  full  array,  the  combination 
of  Events  1 and  2 is  resolved  by  both  processes 
(Figure  38)  ; the  combination  of  Events  1 and  3 is 
resolved  well  by  TSML  but  not  well  by  FKPLOT  (Figure 
39)  because  it  stripped  on  the  wrong  peaks  at  the  right 
side  of  the  figure.  The  combination  of  Events  4 and  5 
fails  again  as  in  the  previous  case  as  can  be  expected 
with  reduced  resolution  (Figure  40).  The  separation 
is  successful  for  the  combination  of  Events  4 and  6 
(F igure  41) . 

Case  2 

One  signal  is  scaled  to  be  three  times  the  ampli- 
tude of  the  other;  no  noise  is  added. 

As  in  the  full  array,  the  combination  of  Events  2 
and  1 (Figure  42)  could  not  be  resolved  well  by  TSML , 
although  there  are  indications  of  the  second  signal 
(dashed  contours).  FKPLOT,  on  the  other  hand,  resolves 
both  signals  with  a little  ambiguity,  a southwesterly 
peak  at  3 dB  below  the  first  peak.  The  combination 
of  Events  3 and  1 (Figure  43)  is  resolved  by  both 
processes.  The  combination  of  Events  4 and  5 
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(Figure  44)  is  resolved  by  FKPLOT  but  not  by  the  TSML 
process.  The  combination  of  Events  4 and  6 (Figure  45) 
is  resolved  by  both. 

Case  3 

Events  combined  with  equal  maximum  amplitudes, 
noise  added.  S/N  ratio  2. 

A brief  summary  of  these  runs  is  as  follows: 
combinations  of  Events  1 and  2 (Figure  46)  and  Events 
4 and  6 (Figure  49)  are  poorly  resolved  by  TSML, 
fairly  well  resolved  by  FKPLOT.  The  other  combinations 
(Figures  47  and  48)  are  not  resolved  by  either  process. 

Case  4 

Events  combined  with  equal  maximum  amplitudes, 
noise  added.  S/N  ratio  1. 

None  of  these  cases  (Figures  50-53)  is  resolved, 
except  the  combination  of  Events  l and  2 (Figure  50) 
by  FKPLOT. 
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F-K  ANALYSES  USING  THE  D RING  AND  CENTER  ELEMENT  ONLY 


General 

This  configuration  has  a very  small  array  diameter 
050  Km)  and  contains  only  five  array  elements  (see 
Figure  72) . This  diameter  is  comparable  to  that  of 
the  seven-element  hexagonal  array  exemplified  by  the 
inner  part  of  ALPA,  which  was  proposed  as  a standard 
for  long-period  arrays.  Although  our  array  has  only 
five  elements,  it  is  fairly  characteristic  of  small 
arrays  consisting  of  only  a few  elements  such  as  the 
seven-element  hexagon.  array.  Because  of  a bad  D2 
trace,  Event  5 and  combinations  involving  it  were 
omitted  from  this  array  analysis. 

Individual  signals  and  noise 

Figures  54  through  58  show  the  f-k  plots  of  the 
individual  events  and  Figure  59  that  of  the  noise. 

The  f-k  spectra  of  the  individual  events  clearly 
demonstrate  the  differences  between  the  appearance  of 
the  outputs  of  both  processes.  The  TSML  outputs  all 
show  a well  localized  peak,  low  or  no  side  lobes, 
while  the  spectral  peak  of  the  FKPLOT  output  is  con- 
siderably broadened,  and  the  side  lobes  are  prominent. 
Although  the  peaks  of  TSML  output  also  broadened  with 
the  reduction  of  the  array  diameter,  this  broadening 
is  relatively  small  compared  to  that  of  the  FKPLOT 
output  (compare  with  Figures  8 through  13). 
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ANALYSIS  OF  COMBINED  EVENTS 


Case  1 

Events  are  equal  in  amplitude;  no  noise  added. 

The  combination  of  Events  1 and  2 (Figure  60)  is 
considered  to  be  resolved  on  both  types  of  outputs, 
although  a special  -4  dB  contour  is  needed  to  show  a 
second  peak  on  the  TSML  output.  However,  the  bulge 
of  the  -6  dB  contour  is  so  conspicuous  that  it 
inevitably  suggests  the  presence  of  a second  signal. 

FKPLOT  resolves  both  signals  with  stripping.  The 
combination  of  Events  1 and  3 (Figure  61)  is  also 
resolved  with  questionable  accuracy.  The  TSML  shows 
a -9  dB  peak  at  the  right  place  for  Event  3.  The 
original  FKPLOT  shows  both  peaks,  but  it  strips  on 
the  wrong  peak  leaving  a very  doubtful  signal  peak. 

The  combination  of  Events  4 and  6 (Figure  62)  is 
resolved  by  FKPLOT  but  not  by  the  TSML. 

Case  2 

One  signal  is  scaled  to  be  three  times  the  ampli- 
tude of  the  other,  no  noise  added  (Figures  65-65). 

This  case  can  be  summarized  as  follows:  the 

combination  of  Events  2 and  1 (Figure  63)  is  not 
clearly  resolved  by  either  method.  Signal  1 on  TSML  out- 
put is  indicated  by  a bulge  in  the  -12  dB  contour.  FKPLOT 
shows  Event  2 clearly,  but  the  stripped  version  is  some- 
what ambiguous.  The  combination  of  Events  3 and  1 


(Figure  64)  is  considered  to  be  resolved  by  both 
processes,  although  an  additional  -7  dR  contour  is 
needed  to  define  a second  peak  on  the  TSML.  The  com- 
bination of  Events  4 and  6 (Figure  65)  is  resolved  by 
FKPLOT  but  not  by  TSML. 

Case  5 

Events  combined  with  equal  maximum  amplitudes, 
noise  added,  S/N  ratio  2 (Figures  66-68). 

None  of  these  combinations  is  resolved  except  for 
Events  1 and  2 (Figure  66)  by  FKPLOT. 

Case  4 

Events  combined  with  equal  maximum  amplitudes, 
noise  added,  S/N  ratio  1 (Figures  69-71). 

None  of  the  combinations  is  resolved. 
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SUMMARY  AND  CONCLUSIONS 


Comparison  of  the  two  segment  maximum  likelihood 
process  (TSML)  and  the  fast  frequency  domain  beaming 
process  ( FKPLOT , FKCOMB)  combined  with  stripping  shows 
that  the  two  processes  are  comparable  in  performance, 
but  the  latter  process  performs  somewhat  better  since 
more  multiple  signal  separations  were  achieved  by 
FKPLOT  as  demonstrated  by  Table  III.  The  advantage  of 
FKPLOT  can  be  attributed  to  the  stripping  process  which 
can  uncover  weak  signals  in  the  presence  of  strong 
signals.  It  is  possible,  however,  to  combine  the  TSML 
with  a procedure  which  is  equivalent  to  stripping, 
such  as  spectral  eigenvector-eigenvalue  analysis,  in 
which  case  TSML  performs  as  well  or  even  better  than 
FKPLOT;  this  is  the  subject  of  a report  now  in  prepara- 
tion. One  clear  advantage  of  the  TSML  process  is  the 
small  dependence  of  the  F-K  spectra  on  the  array 
response.  The  inspection  of  the  TSML  f-k  spectra 
immediately  gives  a clear  idea  about  the  azimuthal 
and  velocity  distribution  of  energy,  while  FKPLOT  gives 
a confusing  picture  full  of  distracting  side  lobes 
which  are  almost  as  large  as  the  main  lobe,  especially 
if  the  array  has  only  a few  widely  spaced  elements 
(see  Figure  64,  for  example).  Another  disadvantage 
of  the  stripping  procedure  is  that  it  can  strip  on  i 
false  signal  consisting  of  the  combined  side  lobes  of 
two  signals  or  sidelobes  combined  by  noise  peaks, 
resulting  in  the  complete  breakdown  of  the  process 


-22- 


/ 


I 


(see  Figures  38  and  63,  for  example).  The  problem 
can  be  minimized  by  stripping  on  peaks  which  fall  into 
the  expected  velocity  range  of  the  signal. 

Another  advantage  of  FKPLOT  is  that  it  seems  to 
be  somewhat  more  stable  in  the  presence  of  noise. 

(This  is  confirmed  by  theory  of  Capon  and  Goodman 
[1970].)  Both  processes  require  a signal  to  noise 
amplitude  ratio  of  about  2 to  perform  adequately. 
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Figure  6.  F;vent  6,  Guatemala,  22  January  197 


igure  7.  Relative  spectral  shapes  for  Events  1 through  6 and  the  noise  sample 


Figure  8.  F-K  spectra  lor  Event 
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Statistical  measures  of  S/.N  ratio  computed  during  the  FKPLOT  runs  are  measures 
of  the  deviation  ot  the  i put  from  a single  plane  wave  rather  than  that  of 

ratios  of  signal  level  to  background  noise. 
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Summary  of  runs  testing  the  relative  resolution 
of  TSML  and  FKPLOT 
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